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ABSTRACT 

Theoretical models of galaxy formation predict that galaxies acquire most of their 
baryons via cold mode accretion. Observations of high-redshift galaxies, while show- 
ing ubiquitous outflows, have so far not revealed convincing traces of the predicted 
cold streams, which has been interpreted as a challenge for the current models. Using 
high-resolution, zoom-in smooth particle hydrodynamics simulations of Lyman break 
galaxy (LBG) halos combined with ionizing radiative transfer, we quantify the cov- 
ering factor of the cold streams at z = 2 — 4. We focus specifically on Lyman limit 
systems (LLSs) and damped Lya absorbers (DLAs), which can be probed by absorp- 
tion spectroscopy using a background galaxy or quasar sightline, and which are closely 
related to low-ionization metal absorbers. We show that the covering factor of these 
systems is relatively small and decreases with time. At z = 2, the covering factor of 
DLAs within the virial radius of the simulated galaxies is ~ 3% (~ 1% within twice 
this projected distance), and arises principally from the galaxy itself. The correspond- 
ing values for LLSs are ~ 10% and ~ 4%. Because of their small covering factor 
compared to the order unity covering fraction expected for galactic winds, the cold 
streams are naturally dominated by outflows in stacked spectra. We conclude that the 
existing observations are consistent with the predictions of cold mode accretion, and 
outline promising kinematic and chemical diagnostics to separate out the signatures 
of galactic accretion and feedback. 

Key words: cosmology: theory - galaxies: formation and evolution, high-redshift - 
radiative transfer 



1 INTRODUCTION 

Recent years have seen the emergence of a new theoreti- 
cal paradigm for galaxy formation, in which galaxies are 
predicted to acquire most of their baryons via cold, fila- 
mentary streams penetrating deep inside dark matter ha- 
los (e.g., Katz et al. 2003; Keres et al. 2005, 2009a). This 
prediction from cosmological hydrodynamical simulations, 
in addition to being supported by idealized ID simulations 
and analytic arguments (Birnboim &; Dekel 2003; |Dekel 
Birnb oim 2006), has now been confirmed by several groups 
using both smooth particle hydrodynamics (SPH) and grid- 
based codes (e.g ., |Ocvirk et aI1|2008| |Brooks et aT1|2009| 
Dekel et al.|2009| ) . Understanding how galaxies get their gas 



ter et al. 2010). Theoretically, the mode of gas accretion 



is a fundamental problem for many reasons. Observationally, 
measurements of the star formation rate density across cos- 
mic time and of the neutral atomic and molecular reservoirs 
show that galaxies must continuously accrete material from 
the ionized intergalactic medium (IGM) in order to sustain 
their star formation ( |Prochaska &; W olfe 2009; Bauermeis- 
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determines the rate at which baryons can settle at the bot- 
tom of dark matter halos, and therefore how rapidly galaxies 
can grow as a function of redshift and host halo properties 
(e.g., |Faucher-Giguere et al.|prep| . Whether baryons flow in 
galaxies in the form of narrow, dense streams or via the cool- 
ing of quasi-spherical atmospheres is also directly relevant 
to the efficiency of feedback processes, and therefore poten- 
tially to the quenching of star formation in massive galax- 
ies and the establishment of the red sequence (e.g., Dekel 
& Birnboim 2006; Keres et al. 2005, 2009b). Finally, the 
cold streams may be connected to several observed phenom- 
ena, including high-redshift clumpy discs (e.g., Elmegreen 
& Eimegreen||2006| |Genzel et al.||2006|), Lya blobs (but 



see Faucher-Giguere et al. 



2010), and high- velocity clouds 



around local galaxies (e.g., Keres &; Hernquist[ 2009). 

In spite of the key role the cold mode is likely to play 
in galaxy formation and evolution, observations have so far 
provided little direct evidence of infalling cold halo gas. 
Thanks to novel techniques, including the use of projected 
quasar-quasar (e.g., Prochaska & Hennawi 2009), quasar- 



galaxy (e.g., Adelberge r et al.|2005 Simcoe et al.|2 006) and 



galaxy-galaxy pairs (e.g., Steidel et al.|2010" ), observers have 
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Figure 1. Projected neutral hydrogen column density (xy projection in Table[l] cm -2 ) around the Bl Milky Way progenitor at z = 2, 3, 
and 4. The dotted contours indicate Lyman limit systems, while the regions corresponding to damped Lya absorbers are shown by the 
solid contours. The point deepest in the potential is at the center of each panel and the dashed circles mark the virial radii. A gas 
smoothing length of 27 proper pc is achieved at z = 2, making this one of the highest-resolution simulations of such a halo to date. 



recently made impressive progress in probing the halo gas 
content of high-redshift galaxies. While the spectroscopic ab- 
sorption measurements show strong evidence for ubiquitous 
galactic outflows at z ~ 2 — 3, they do not show clear sig- 
natures of the predicted accretion streams, which has been 
interpreted as a challenge for current models (e.g., [Steidel 



et al.|[20T0 ). Exactly what these observations tell us about 
the cold mode has however been ambiguous, as we are only 
beginning to quantify how the predictions for the accretion 
dynamics should manifest themselves in observables. The 
opportunity that is presented to us by this state of affairs 
can hardly be understated, as it is a unique case of a nearly 
ab initio physical prediction in galaxy formation that has 
yet to be tested, but for which observations are rapidly ap- 
proaching truly discriminating power. It therefore behooves 
us to robustly quantify the predictions of our models. 



In|Faucher-Giguere et al." (2010), we made a step in this 



direction by calculating the Lya cooling emission from cold 
streams, showing that the treatment of self-shielded gas in- 
troduced orders-of-magnitude uncertainties in most previous 
work, and concluding that pure accretion cooling is unlikely 
to dominate the Lya emission of the bright Lya blobs now 



routinely imaged (e.g., Matsuda et al. 2010), although it 



could be important for some fainter sources (e.g., Rauch et 
|al.|2008| . In this Letter we address absorption spectroscopy 
as a different observational probe of the cold streams. Us- 
ing high-resolution cosmological simulations of a Milky Way 
progenitor and of a typical LBG-mass halo, we quantify the 
covering factor of the cold streams in Lyman limit (LLSs) 
and damped Lya systems (DLAs). We show that the cover- 
ing fraction of such dense gas decreases with time and that it 
is sufficiently small at z ~ 2 that the presence of cold streams 
in the halos of Lyman break-selected galaxies is consistent 
with current observations, consisting principally of stacked 
spectra, that are dominated by outflow signatures. 



2 SIMULATIONS AND RADIATIVE 
TRANSFER 

2.1 Hydrodynamics 

Our simulations use a modified version of the GADGET 
2 cosmological code (iSpringel 2005). The gas dynamics is 



calculated using an SPH algorithm that conserves both en- 
ergy and entropy (Springel &; Hernquist|200"2 ). The modifi- 
cations with respect to the public version of the code include 
the treatment of cooling, the effects an uniform ultra-violet 
background (UVB), and a multiphase star formation algo- 
rithm as in Springel & Hernquist (2003). The thermal and 
ionization properties of the gas are calculated including all 
relevant processes in a plasma with primordial abundances 
of hydrogen and helium following Katz et al. ( 1996 ), using an 
update of the [Haardt &; Mada"u| ( |1996| UV background in- 
cluding galaxies and quasars. Self-shielding effects are mod- 
eled in post-processing, as described in the next section. To 
achieve ultra-high resolution, we 'zoom in' on individual ha- 
los within a (10 /h comoving Mpc) 3 simulation box and fol- 
low the local gas dynamics at a resolution refined 64 x in 
mass. 

Our main simulation of a Milky Way progenitor (labeled 
Bl) is identical to the one analyzed by |Keres &; Hernquist 
(2009), but with a mass resolution higher by a factor of 
3.375 because we used 192 3 dark matter particles in the un- 
refined volume instead of 128 3 . The minimum achieved gas 
smoothing length is 27 proper pc at z — 2, with a Plummer- 
equivalent gravitational force softening of 92 [(1 + z)/3]~ 1 
proper pc. The halo mass ~ 3 x 10 11 Mq at z = 2 cor- 
responds to a low-mass LBG, and the simulation is among 
the highest-resolution cosmological SPH realizations of such 
a halo to date. The dark matter particle mass of this simu- 
lation is 2 x 10 5 Mq, and gas particle mass is 4 x 10 4 Mq. 
We also study a zoom-in simulation of a more typical LBG 
halo (labeled Al) of mass ^ 9 x 10 11 M at this redshift, 
with a spatial resolution lower by a factor of 1.5. We as- 
sume a flat ACDM cosmology with Q m = 0.27, ^ b = 0.044, 
h = 0.7, as = 0.8, and n s = 0.95, consistent with the latest 
WMAP analysis ( |Komatsu et al.|[2010| . Galactic winds are 
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Small Cold Streams Covering Factor 

Table 1. Covering Factors of Cold Mode Lyman Limit Systems and Damped Lyo: Absorbers Around the Bl Milky Progenitor 



3 



z 


M h 


-Rvir 


LLS(< 0.5i? vir ) 


DLA(< 0.5i? vir ) 


LLS(< i? vir ) 


DLA(< R vir ) 


LLS(< 2i? vir ) 


DLA(< 2i? vir ) 




M© 


prop, kpc 


% (xy/xz/yz) 


% (xy/xz/yz) 


% (xy/xz/yz) 


% (xy/xz/yz) 


% (xy/xz/yz) 


% (xy/xz/yz) 


2 


3.2 x 10 11 


74 


32 / 40 / 23 


10/9/10 


10 / 15 / 11 


3/4/3 


4/4/4 


1/1/2 


3 


2.5 x 10 11 


52 


33 / 23 / 45 


17 / 34 / 17 


16 / 16 / 23 


6/12/5 


8/5/10 


3/3/2 


4 


7.4 x 10 10 


27 


50 / 29 / 51 


25 / 53 / 25 


32 / 27 / 32 


11 / 25 / 10 


17 / 15 / 14 


6/13/6 



The three covering factors given for each case correspond to projections along different Cartesian axes (on the xy, xz, and yz planes). 



neglected in this work to isolate the accretion streams, but 
will be studied in future work. We define the virial radius, 
i?vir, as the radius enclosing a mean overdensity of 180 times 
the mean matter density. 



2.2 Radiative Transfer 

The covering factor of the cold streams can only be meaning- 
fully quantified for specific boundary criteria. For absorption 
statistics a natural definition is tied to the column density of 
the material. We limit our attention here to hydrogen, and in 
particular to the dense LLSs (10 17 2 ^ Nm ^ 2 x 10 20 cm" 2 ) 
and DLAs (iVm > 2 x 10 20 cm -2 ) that are often associated 
with proto-galaxies and their outskirts. By definition, these 
systems can shield themselves from the ionizing radiation of 
the cosmic background and local sources. To model the ion- 
ization balance, we use the ray tracing code introduced in 



Faucher-Giguere et al. (2010). The UV background is mod- 
eled by following 512 rays originating from the 6 faces of 
a cubic volume of side length 1/h comoving Mpc centered 
on the halo of interest. We have improved this code for the 
present work by including frequency-dependent attenuation 
of the radiation field, and by allowing for the possibility of 
enhanced ionizing flux from local sources by radially casting 
rays outward. 

Assuming that the central galaxy emits ionizing pho- 
tons in proportion to its star formation rate, for a Salpeter 
IMF and for an escape fraction of 5%, we find that ionization 
from local sources has a very small impact on the covering 
factor of surrounding LLSs, and is completely negligible for 
DLAs. This is easy to understand, as the optical depth of 
DLAs to ionizing photons exceeds 10 3 , so that the bulk of 
their gas does not see changes in ionizing luminosity unless 
they are extreme. Because the covering factors of LLSs and 
DLAs are insensitive to local ionization, we chose not to in- 
clude it in the results presented here, which has the benefit 
of ensuring that we are providing upper bounds (in the ab- 
sence of additional physics), and therefore strengthens our 
main point that they are small for the cold streams, espe- 
cially at z rsj 2. To avoid artificial collisional ionization of 
the interstellar material (ISM) due to the effective temper- 
atures carried by the multiphase particles, we assume that 
these are fully neutral, up to a maximum nm — 1 cm -3 to 
model the conversion to H 2 (e.g., |Schaye| 2001 ). 



3 RESULTS 

Figure [I] shows the projected HI column density distribution 
around the simulated Bl Milky Way progenitor at z — 2, 3, 
and 4. Table [l] summarizes the properties of the system at 
these redshifts and quantifies the covering factors of LLSs 
and DLAs. At each redshift, we provide the covering factors 
averaged within 0.5, 1, and 2 virial radii, and for three mutu- 
ally orthogonal projections labeled according to the Carte- 
sian projection plane (xy, xz, or yz), where the xy projection 
is the one shown in Figure [I] The covering factors are eval- 
uated exactly by considering all the pixels of the projected 
radiative transfer grids. Providing the covering factors nor- 
malized to fractions of the virial radius gives a simple way 
to scale the results to halos of different masses. However, 
the availability of background sources in observations is in- 
dependent of the virial radius of the foreground galaxy. In a 
forthcoming study, we will also provide numbers as a func- 
tion of absolute distance for samples representative of ob- 
servations, but it is already clear from Figure [I] that much 
of the covering factor evolution within fixed fractions of the 
virial radius is driven by the evolution of the virial radius 
itself, and that the covering factors within a fixed distance 
do not evolve as strongly. For reference, the covering factors 
of DLAs [LLSs] between fixed radii of 10 and 100 proper 
kpc, for the xy projection of the Bl halo, are (1, 2, 3)% 
[(7, 8, 10)%] at z= (2, 3, 4). 

At z = 4, where the streams are prominent and the 
virial radius (27 proper kpc) is compact, 10 — 25% of the area 
within i? v ir is covered by a DLA, depending on the sightline. 
By z = 2, however, this fraction is reduced to 3 — 4% and 
even so most of it arises from galactic material rather than 
from the accretion streams, which are seen to actively frag- 
ment (excising the inner 10 proper kpc in radius, the DLA 
covering factor drops to 1 — 3%, depending on the sightline). 
Within 2 i? v ir, the DLA covering factor plummets below 1% 
even including the galaxy. The filament fragmentation in 
this high-resolution simulation of a ~ 3 x 10 11 M© halo at 
z = 2 is consistent with previous, lower-resolution work that 
showed that streams survive mostly in halos below the tran- 
sition mass ~ 2 — 3 x 10 11 M©, or at z > 2. At lower redshift, 
cooling and possibly Ray leigh- Taylor instabilities start op- 
erating in the infalling gas, which is increasingly surrounded 
by a hot medium, and seed the formation of cool clouds as 
in |Keres fc Hernquist] ( |2009| ). 

Figure [2] shows the HI column density map around the 
Al halo, of mass « 9 x 10 11 M at z = 2. While the resolu- 
tion of this simulation is poorer by 50% spatially, the halo 
mass is exactly the average mass probed by the LBG sample 
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z=2 



similar to HI (e.g., Steidel|1990 




lei 1 



Figure 2. Same as Figure[T] but for the Al halo of mass 9 x 10 11 
Mq at z = 2, corresponding to an average-mass LBG in the 
|Steidel et al] ( |2010| > sample. 



oflSteidel et al.l (120101) at the mean redshift (z) = 2.2, and 
allows us to test how the results above scale to this mass. For 
the random projection shown, the covering factor of DLAs 
within i?vir at z = 2 is 4%, i.e. the same as for the xz projec- 
tion of the Bl halo in Figure [I] and again arises mostly from 
material close to the galaxy. The LLS covering factor of 11% 
within i?vir is also similar to the Bl case. This result is con- 
sistent with the study of Hong et al. (2010), who find that 



the DLA covering factor within R v i r in simulations without 
outflows is relatively constant over a fairly wide mass inter- 
val in the LBG range. The small covering factors of accreting 
material we find are therefore generic predictions for LBGs 
at this redshift. In simulations with outflows, |Hong et al.| 
(2010) show that the DLA covering factor increases in high- 
mass halos, but their kinematics suggest that most of the 
increase owes to absorption by wind ejecta. Our point re- 
garding the accretion streams is therefore unaffected to first 
order, although it is conceivable that interactions with out- 
flows might affect their dynamics somewhat. Between the 
fixed radii of 10 and 100 proper kpc, the covering factors of 
DLAs and LLSs are almost identical, at 4% and 12% respec- 
tively, but higher than the corresponding numbers in the Bl 
case owing to the overall difference in halo size. 



4 DISCUSSION 

4.1 Comparison with Observations 

A key question is whether the existing observations are con- 
sistent with cold accretion. While we have quantified the 
covering factors of LLSs and DLAs around z = 2—4 galaxies, 
direct observational measurements of these quantities are 
not yet available in most cases. Our comparison is therefore 
somewhat qualitative at this time, as more accurate infer- 
ences will require larger simulation samples and the explicit 
modeling of metal absorbers. LLSs and DLAs are however 
closely connected to low-ionization metal species such as OI, 
Sill, CII, Fell, and Mgll, which have ionization potentials 



Wolfe et al.|2 005). It is also 
clear from from Figures [1] and 2] that the column densities 
of neutral gas drop sharply outside the self-shielded regions. 

The most exhaustive constraints on the circum-galactic 
medium of z ~ 2 — 3 galaxies were recently reported by 
Steidel et al. ( 20101) , based on spectroscopic observations of 



a large number of LBGs, including 512 close angular galaxy- 
galaxy pairs. While both the blueshifted interstellar absorp- 
tion and the redshifted Lya emission in direct sightlines 
provide strong evidence for high- velocity outflows, these au- 
thors concluded that there is at present little evidence of 
infalling cool material. In fact, they showed that a simple 
outflow model calibrated from direct sightlines also repro- 
duces the observed transverse statistics in galaxy pairs. Un- 
like background quasars, the vast majority of LBGs are too 
faint to yield spectra of sufficient resolution and signal-to- 
noise ratio to study individual absorbers. The main results of 



Steidel et al. (2010) are consequently based on stacked mea- 



surements. The authors also argue that the principal metal 
lines they consider are strongly saturated, so that the ob- 
served equivalent widths are good proxies for their covering 
factors and insensitive to metallicity. If that is the case, then 
infalling material would leave negligible traces in composite 
spectra of transverse absorption whenever its covering fac- 
tor is much smaller than that of outflows. Since the mean 
redshift of Steidel et al.'s foreground galaxies is (z) = 2.2, 
this hypothesis is well supported by the small cold stream 
covering factors we predict at z ~ 2. As shown by Hong et al.| 
( |2010| ), outflows can in contrast fill a substantial fraction of 
the virial radius with self-shielding systems. Note that we do 
not confirm the prediction of a larger, ~ 25% covering factor 
of DLAs between r adii 20 and 100 proper kpc predicted by 
Dekel et al.] ( |2009| ) for 10 12 M halos at z = 2.5. The sim- 
ulation analyzed by those authors had a spatial resolution 
than lower than ours by a factor ~ 15 at this redshift, and 
they did not perform a radiative transfer calculation but in- 
stead included all the gas with temperature T < 10 5 K in 
their column densities. Both effects go in the direction of 
overestimating the covering factors. 

The above argument does not apply to direct single 
sightlines, in which velocity offsets v\ os ^ might be more 
directly associated with infall and in which absorption is 
most likely to occur close to the center. However, |Kimm| 



et al. ( |2010| recently demonstrated that the infalling mate- 



rial (with velocities of order the halo circular velocity, v c 
is in general completely masked by absorption by the ro- 
tating ISM of the galaxy itself. With regards to Lya emis- 
sion, ^^^^^Gi^^Te^^]^010^ showed that accretion via 
cold streams (again because of their small covering factor, 
which allows the photons to take paths of least resistance be- 
tween them) can only produce a minor blue enhancement. 
Winds with larger covering factors therefore drive the radia- 
tive transport and yield redshifted lines. 

We stress that the cold accretion paradigm is consistent 
with the presence of galactic winds. We in fact fully expect 
the real high-redshift galaxies to drive powerful winds, as 
appears required to explain a number of observations be- 
yond those discussed above, including the low-end of the 
galaxy luminosity function and the pollution of the IGM 
with metals. The omission of galactic winds in the simula- 
tions analyzed here simply reflects the fact that, whereas 
simulating accretion from the cosmic web is a well-defined 
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problem, we do not yet know how to robustly model winds. 
We have therefore deliberately excluded them to isolate the 
signatures of the cold streams. This is a reasonable ap- 
proximation, as the small covering factor of the dense cold 
streams implies that it is difficult for the more dilute winds 
to strongly affect them (e.g., Faucher-Gigu ere et al.||prep ). 



4.2 Inflows vs. Outflows Diagnostics 

The present lack of direct observational evidence for cold 
accretion and the small predicted stream covering factors, 
especially relative to the ubiquitous galactic winds at high 
redshift, highlight the challenge of testing the theory. We 
conclude by outlining ways in which future studies might be 
able to reveal the cold accreting gas. 

The most direct diagnostics of accretion versus winds 
are the kinematics and chemical composition of cool ab- 
sorbers. In transverse spectra, the direction of motion of 
the absorbers is usually ambiguous because their phys- 
ical positions are unknown. However, the magnitude of 
the velocity offset provides a discriminant, as material 
infalling from the IGM should have a velocity v c i rc ~ 



245 km s 



-1 ( M h \ 1/3 ( l+z\ 
\10 12 M J \ 3 ) 



,1/2 



whereas outflows are 



io 12 m / 

observed to have velocities up to 800 km s _1 or more (e.g., 
Steidel et al.|2 010). A caveat is that wind ejecta can in prin- 
ciple fall back onto galaxies with velocities ~ v c i rc , if their 
launch velocity is modest or if the galactic potential well is 
sufficiently deep (e.g., Oppenh eimer et aT||2010 ). Metallic- 
ity may then be the only means to convincingly associate a 
particular absorber with intergalactic accretion. By carefully 
modeling the absorbers in their high signal-to-noise, high- 
resolution quasar spectra, Simcoe et al. ( 2006 ) were able to 



derive the metallicity of individual systems near galaxies, 
and link them to feedback processes based on their near so- 
lar abundances. By contrast, the IGM metallicity is ~ 1/700 
of solar at z = 2.5 (e..g, Simcoe et al.|2004 ), although mate- 
rial infalling onto galaxies could be substantially more pre- 
enriched owing to dwarf galaxies in filaments, and it is an 
important question for future modeling exactly what are the 
typical and maximum metallicities expected for cold mode 
accretion. 

A practical issue is that most of the information cur- 
rently available on absorbers within the halos of LBGs is 
from close galaxy-galaxy pairs (Stei deTet al.|2010 ). Because 
galaxies are faint, it is in general not possible to model indi- 
vidual absorbers and infer their metallicity. The holy grail 
would be a large sample of close quasar-galaxy pairs from 
which to build up the kinematic and metallicity statistics 
of halo absorbers, but the rarity of quasars limits the use- 
fulness of this approach. Deep spectroscopy of gravitation- 
ally lensed galaxies would also be very helpful, and next 
generation 30-m class telescopes will greatly improve the 
resolution and signal-to-noise ratio of galaxy spectra, and 
therefore enhance the information accessible via the galaxy- 
galaxy technique. Sufficiently sensitive integral field spectro- 
graphs would be ideal instruments to map the extended gas 
distribution around LBGs. 

These observational proposals are however challenging, 
and underscore the critical role to be played by more de- 
tailed theoretical modeling. In analogy to how statistics from 
a large number of galaxy-galaxy pairs partially make up for 



the rarity of suitable quasar-galaxy pairs, more statistical 
predictions will connect our models better with the exist- 
ing data. For example, it should be possible to theoretically 
compute the average equivalent width of different ions versus 



impact parameter, as measured by | Steidel et al. (2010) from 
galaxy-galaxy pairs, and use this information to constrain 
the models. For this exercise to be fruitful, it will however 
be necessary to improve our modeling to consistently in- 
clude the galactic winds that produce most of the measured 
absorption, the attendant chemical enrichment, and to sim- 
ulate a sizable sample of halos at sufficient resolution. We 
plan to address these problems in future work, as well as 
to investigate the constraints that can be derived on high- 
redshift outflows ( Faucher-Giguere et al.Jprep ). In the mean 
time, it is becoming clear that the signatures of cold accre- 
tion are quite subtle, and that accurate predictions will be 
needed to make sense of the data. 
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